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Abstract 

Correlators of the octet baryons in the hot pion gas are studied in the framework of the 
QCD sum rule. The condensates appearing in the OPE side of the correlators become T- 
dependent through the interaction with thermal pions. We present an explicit demonstration 
that the 0(T^)-dependence of the condensates is completely compensated by the change of 
the pole residue and the tt + B ^ B' scattering effect in the spectral functions. Therefore 
the baryon masses are constant to this order, although {uu)t — {uu)o{l — which is 

consistent with the chiral symmetry constraint by Leutwyler and Smilga. 



1 Introduction 



Modification of tlie liadronic properties at finite temperature has been receiving greater 
attention in connection with the relativistic heavy ion colhsions such as RHIC and LHC 
Although the finite-T properties of various mesons have been extensively studied in the 
literature using QCD Sum Rules (QSR) [0, H, and effective theories for QCD 0, H, ^ [^l, 
baryon properties have not been seriously investigated except for the nucleon p, ^ |T0[ . In 
principle, the change of baryon masses should manifest itself in the yield modification and 
the threshold behaviour of lepton pairs coming from the baryon-antibaryon annihilation in 
the relativistic heavy ion collisions. Therefore it is a pressing issue to have a QCD prediction 
for the finite-T behavior of the baryon correlators. 



Since the work by Shifman, Vainstein and Zakharov |TT| , the QCD sum rule method has 
been extensively used as a systematic tool to study various resonance properties based on 
QCD (see [0 for a review). Especially, the masses of octet and decuplet baryons have 
been well reproduced in terms of the vacuum condensates [1T3| , ^ 0. Hence it is worth 
while to extend this strongly QCD motivated phenomenology to study finite-T behavior of 
baryon correlators. In this paper we shall analyze the correlators of the octet baryons (N, 
A, S and S) at finite temperature. The source currents we will investigate are the following 
interpolating fields which have been used in QSR: 



^^(x) = e'^''u^{x)C-fy{xh,^^d%x), (1.1) 

(1.2) 



V'^ix) = y^e'^^'=[M'^(x)C7;.s^(x)757M^'=(a^)-^^"(a:)C7^s''(x)757X(a^) 

r^^(x) = e''''u''{x)C^y{x)^,^,s'ix), (1.3) 
r^^(x) = -e^''s''{x)Cj,s\x)j5lf.u'{x), (1.4) 

where a, b, c are the color indices and C is the charge conjugation matrix. Although one 
can choose other sets of source currents for the octet baryons, it has been known that the 
above combination gives the best description of the octet baryon masses. 

A useful quantity for studying the finite-T behavior of baryons is the retarded correlation 
function for the above fields [|l7l : 

Ii%{uj,q,T) = i j d^xe''^^e{x,){Ux)np{Q) + np{Q)r^a{x))T^ (1.5) 

where q = {lo, q) is the external four momentum, {O) = Tr {Oe~^^'^^ /Z {Z = Tr(e~^/^)) 
is the thermal average for the operator O with H being the QCD hamiltonian, and a and f3 
are the spinor indices of the interpolating fields ( p..l|) -( p!^ . The retarded correlator satisfies 
the dispersion relation with an appropriate spectral function p{u!, q,T): 

^ J-oo u — 00 — te 

In the QCD sum rule method, one applies the operator product expansion (OPE) to the left 
hand side of the correlator in the deep Euclidean region = — oo and adjusts the 
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resonance parameters (resonance masses, pole residues, width of the resonances, continuum 
threshold, etc) in the spectral function p so as to reproduce the OPE side of the correlator. 
This procedure provides us with an expression for the resonance parameters in terms of the 
condensates. Therefore, in the QSR approach, the change of the condensates at finite-T 
naturally causes a response of the baryons as a change of their properties at finite-T. For 
example, one might naively expect that the nucleon mass will drop at finite temperature as 



the chiral order parameter decreases . 

On the other hand, Leutwyler and Smilga ||^ considered the nucleon correlator in a ther- 
mal pion gas and showed that the nucleon mass does not have 0(T^)-dependence: In the 
thermal pion gas, the T-dependence of the nucleon correlator is associated with the pion- 
nucleon forward scattering amplitude. It becomes zero in the soft pion limit (Adler's zero) 
and thus the real part of the self-energy becomes zero at the pole position of T = 0. The 
same statement is also true for all the octet baryons. 

The purpose of this paper is to demonstrate how this seeming contradiction can be 
reconciled in the framework of the finite-T QCD sum rules. We organize the QSR at T 7^ 
for the nucleon in the dilute pion gas . This approximation should be valid at relatively 
low temperature (T < 150 MeV) well inside the confined phase. In this framework, the 
condensates appearing in the OPE side of the correlators receive T-dependence through the 
pion matrix elements of the same operators. We pay particular attention to the effect of 
the TT -|- iV — > iV scattering in the phenomenological side of the QSR as was suggested by 
Eletsky [^]. Then it can be shown that the 0(T^)-dependence of the condensates is exactly 
compensated by this scattering and the mass does not have 0(T^)-dependence as was found 
by Leutwyler and Smilga. Thus the naive expectation motivated by the loffe's mass formula 

1 /3 

for the nucleon, Mjv(T) ~ (— 2(27r)^(-uu)r) , does not work at finite-T. The situation is 
analogous for the other octet baryons. 

The construction of this paper is the following: In section 2, after briefly summarizing 
our finite-T QSR for the case of the nucleon along the line of 0, we present an explicit 
demonstration that the nucleon mass does not have a 0(T^)-dependence in the QSR. The 
discussion for the hyperons is similar. So we will not repeat it, but present some of the 
formulas in the Appendix. In section 3, we shall analyze the baryon correlators in the pion 
gas in terms of PCAC without using QCD sum rules. We shall see that the octet baryon 
correlators at finite-T can be written in terms of the correlators at T = with T-dependent 
coefficients and our OPE expression derived in section 2 is consistent with those relations. 
This explains why the pole positions of the octet baryon correlators do not have 0(T^)- 
dependence. Section 4 is devoted to a brief summary and outlook. 



2 QCD Sum Rules for Nucleon at T 7^ 

2.1 OPE at T 7^ 

We shall first discuss the QCD side of the correlator. For simplicity, we study the nucleon 
at rest, i.e. q = Q. In the region, Q"^ = —q^ = —uP' > 0, the retarded correlator ( |1.5|) is 
identical to the causal correlator: 

Uo.p{u;,q,T) = i f d'x e'^'^T {r]^{x)r]p{0)))T. (2.1) 
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At T 7^ 0, Ilai3{q,T) can be decomposed into the three scalar components: 

Il^^iq,T)=U,iq,T)5^p + U2iq,T)^^f, + Mq,T)ii^f3, (2.2) 

where is the average four-flow velocity of the medium equal to = (1, 0, 0, 0) in the rest 
frame. In the deep Euchdean region oo, one can apply OPE to n(g,T), which, for 112 

as an example, can be schematically written as, 

U,{q,T) = T.Cr''{q,f^)Oi^:Uf') 

i 



+ ^ ((0'"°'>T + '-^{O'^fh + ^f^io'^lh) + ■ • (2.3) 

where the i-th local operator 0^^^:!^j^^{fj,) renormalized at the scale /i has dimension-d and 
spin-s with s Lorentz indices and Cf ^ (g, /i) is the corresponding Wilson coefficient. As a 
complete set of the local operators in { \^.'S\ ), one can always choose symmetric and traceless 
operators with respect to all Lorentz indices. We will hereafter assume this symmetry con- 
dition for all nonscalar operators. In the above equation the following two features peculiar 
to T 7^ are implemented : 

(i) T-dependence of the correlators appears only as a thermal average of the local op- 
erators in the OPE as a consequence of the QCD factorization. This is indeed natural if 
we note that such a soft effect should be ascribed to the condensates (C*)t as long as the 
temperature is low enough compared to the separation scale /x, i.e. T <^ n <^ Q. 

(ii) At T 7^ 0, there is no Lorentz invariance due to the presence of the thermal factor 
Q-H/T ^ and hence nonscalar operators survive as condensates: 

{0'^.,■■■^.X^'))T = (mmi ■ ■ -M;,, - traces) a' (/i,T). (2.4) 

At relatively low temperature in the confined phase, the system can be regarded as a non- 
interacting gas of Goldstone bosons (pions). In this approximation, T-dependence of the 
condensates can be written as 

{0^)t ^ {O^) + j:J ^^,{7r''{pW\7r^{p))nB{e/T), (2.5) 



where e = ^Jp^ + m^, a denotes the isospin index, ub^x) = [e^ — 1] ^ is the Bose-Einstein 
distribution function, and (■) is the usual vacuum average. Here we have used the covariant 
normalization for the pion state: (vr°(p)|7r*(p')) = 2e{27r)^6"'^6^{p — p')- Thus we need pion 
matrix elements of the local operators appearing in the OPE to carry out the finite-T sum 
rule. For the scalar operators C)^'^'^), we can apply the soft pion theorem 



{n^ip)\0\Ap)) = ^(0| ^5,0 |0) + ( -p^ ) , (2.6) 



■1 



m. 
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^HAD, 
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where Ah ad is a typical hadronic scale of O ( 1 Ge V) and JF^ is the isovector axial charge 
defined by 

^5^ = J rf'xg(x)7o75Yg(a:). (2.7) 

Pion matrix elements of the nonscalar operators are associated with the pion structure func- 
tions measurable in the hard processes (deep inelastic scattering, Drell-Yan, direct photon 
production etc.). Experimentally, however, only twist-2 part of these matrix elements are 
known to some extent, and therefore we have to await future precise measurements of the 
twist-4 pion structure functions to carry out the satisfactory QSR at finite-T. Pion matrix 
elements of these operators read 

{Hp)K,...Jn{p)) ~ {p^, traces)A,(/i), (2.8) 

and therefore the contribution of nonscalar condensates becomes the effect of O(T^) or higher 
in T (m^r ~ T is assumed) as is easily seen by inserting ( p.8|) into (p^ ). In [|| we found 
that these effects can be neglected below T = 160 MeV. We shall consistently ignore these 
condensates in this work. For the same reason, n3(g, T) also becomes higher order with 
respect to temperature. The effect of the heavier resonances (K, t] etc) was also found to be 
negligible at T < 160 MeV because of the suppression coming from the distribution function 

Applying the soft pion theorem to scalar operators appearing in the OPE of the nucleon 
correlators, we get 

nf(«,r) = -L(s„,) A - ,2,„(q2)_ (2.9) 

nJ'te.T) ^ ^,W)-5^eG^)>nWV?|p2. (2.10) 

where C = {T^ f^)Bi{m^/T) with 

In ( p.9| ) and ( p.lQ ) we note the following points: 

(i) Although we discarded the terms proportional to m„ and because of their smallness, 
we included the pion mass correction coming from the Bose-Einstein distribution, since it 
appears in the form of m^/T in Bi. Bi approaches 1 at -C T, while it is strongly 
suppressed as Bi ~ q-'^^I'^ at ^ T. 

(ii) The chiral condensate {uu)t changes as {uu)o{l — C/8)- Its T-dependence was cal- 
culated by the chiral perturbation theory including up to O(T^) effects |jT^, which gives the 



same coefficient as above for the O(T^) effect. In the temperature range T < 160 MeV, both 
calculations agree well within 5 %. 

(ni) By using the QCD trace anomaly, T-dependence of (^G'^)t can be estimated 0. 
It gives a negligible change of the condensate at finite T (0.5% at T = 200 MeV). We thus 
ignored its T-dependence in (|2.10|) . 
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(iv) The four-quark condensate in the nucleon channel turned out to be T- independent, 
which is quite different from the behavior of the square of the chiral order parameter. Here 
we used the vacuum saturation assumption, i.e., ((-ur-u)^) {uu)"^ after applying the soft 
pion theorem, as is usually adopted in the vacuum QSR. The calculation of the four-quark 
condensates is somewhat tedious, so we shall present a nonfactorized form of the four-quark 
condensates for the octet baryons in the Appendix. There one sees that the T-dependence 
of the four-quark condensates is different in different channels. For example, the four-quark 
condensate in the nucleon channel is T- independent, while it is proportional to 1 — C/6 in 
the III structure in the E-channel (See also ( |A.11| )- (|A.15|) in the Appendix). 



2.2 Spectral function 

In the pion gas, the spectral function for the nucleon current acquires a contribution from 
the scattering with the thermal pions: tt + N ^ N , A etc. The effect of these scatterings 
has to be taken into account in the spectral function before the change of the condensates is 
ascribed to the shift of the pole position. For the purpose of identifying these contributions, 
we introduce the expression for the spectral function [p!7[] : 

Pap{uJ,q,T) = -Imnf^(u;,g,r) 

TT 

= ^(27r)3^(n|r/,(0)|m)(m|r7^(0)|n) 

n,m 

X (e--/^ + e-^-'^)5{u - Urnn)S^'\q - p_), (2.12) 

where the states |m) and \n) have the four momentum {em,Pm) ^'^^ (^n,P„), respectively, 
and Umn = £m - gn, Pmn = Pm - Pn- If ^6 put \n) = |0), |m) = \N{p)) ( |m) = |0), 
\n) = \N{p))) in ( |2.12| ), this is the contribution from the nucleon (anti-nucleon) at T = 0. 



By introducing the nucleon and the anti-nucleon spinor by the relation 

{0\r]^iO)\N{p)) = XNU^ip), {0\UO)\N{p)) = X^Vaip) (2.13) 

with the normalization u{p)u{p) = 2Mn and v{p)v{p) = —2Mn, we get for this contribution 

piu, q) = %4 + Mm) {6{uj - p') - 6{u; + /)) (2.14) 



2po 

with p° = \J q^ + M^. Using the Borel sum rule method, we can study the T-dependence of 
the mass Mj^{T), of the pole residue X^iT), and of the continuum threshold So(T). However, 
it is difficult to incorporate the effects of the width and the scattering contribution induced 
in the pion gas. Thus we shall first list up the effects which should not be associated with the 
change of the above three resonance parameters. Then we put these additional structures 
at T 7^ in the spectral function when we carry out the Borel sum rule. 

(i) TT + N ^ N {tt + N ^ N) contribution; \n) = \iT{k)), \m) = \N{p)) (|n) = \N{p)), 
\m) = |7r(A;))): By taking into account the Bose symmetrization among pions which equally 
fill both \n) and |m), but do not interact with the nucleon current, we arrive at 

p'^+^-^(g,r) = (2vr)3/-4^n^(fc7T) ^ 



(27r)32A;0 ' 'J {2Tif2p^ 
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X 



spin=±l/2 a 



- p° + {q-p + k) + 6{u-p^ - A;°)(5(=^) {q-p-k 



+ E E(^WI^«(o)k'^(fc))K(^)l%(o)|iv(p)) 

spin=±l/2 a 

X (5{uj + / + (q + p + fc) + + p° - U^)5^^^ {q + p-k, 

(2.15) 

Here we discarded the thermal factor for the nucleon (Fermi distribution function), since 
X/(^gAfiv/T _ 1) ~ at T < 200 MeV. We have to include two kinds of contribution to the 
matrix element (7r"(A;)|77„(0)|A^(p)) ]19[. 

(a) Direct coupling of r]^ to the pion (Fig. 1): This contribution can be calculated by 
applying the soft pion theorem: 



{^"mvmiNiP)) = ^(o|[-^5^^s(o)]|iv(p)), 

In 



[^5,^(0)] 



22 



75^? 



12 



(2.16) 



(2.17) 



where r]^ is the proton current defined in ( |1.1[ ) and r/" = —e°'^'^d"'{x)C'y^(f{x)'y5'y^u'^{x) is the 
neutron current. Using (|2.16| ) and ( |2.1?| ) in (p.l5| ) and putting q = 0, we arrive at Q 



-,nB{k'/T) X 



(27r)32P2p' 

75(2^ + Mn)i^ (5{uj - / + fcO) + 5{uj - / - e)) 
+ 75(2^ - M;v)75 ('^(c^ + / + A;°) + 5(u; + p° - A;°) 



(2.18) 

Since we used the soft pion theorem in calculating the matrix element in ( ^.l(j| ), it suffices 
to consider the soft pion limit k = (fc°, fc) — > in (|2.18| ). Then ( |2.18| ) becomes 



32M 



N 



{b{uj - Mm) - 5{uj + Mn)) {4 - Mn) . 



(2.19) 



Equation ( |2.18| ) was derived by Eletsky |]^] taking the imaginary part of the corresponding 
retarded correlator.^ There it was shown that (|2.18| ) has a localized structure around the 
nucleon and anti- nucleon poles, which can be well approximated by ( p. 191) . 

^Here p has k or ~k corresponding to two (^'^■^^-functions in ( |2.15| ). 

^The result given in eq. (11) of |^ is two times larger than the one given in (2.18) by mistake. I thank 
V. L. Eletsky for correspondence to clarify this point. 
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(b) Coupling of rj^ to the nucleoli which interacts with vr (Fig. 2): This contribution can 
be calculated by using the vertex 



(vr(A:)|r^(0)|iV(p)) = A^(7r(fc)|^^(0)|iV(p)) = p , / Mp), (2.20) 

where we assumed the ir — N — N interaction lagrangian as 

= ^N^^l.rNd^Ti. (2.21) 



Inserting ( |2.2(j| ) into ( p^.lSj ), one obtains in the g = limit: 



X [{4 + Mm) ^75 (^ + Mn) ^75 {4 + M^) 

X (5(cu - / + A;°) + 5(cu - / - k^)) 

+ (-^ + Mat) 1^75 - Mm) (-4 + Mm) 

X {6{uj + + k°) + 6{uj + p° - k'^))] (2.22) 



^?a'A^C \ ^^^^ _ _ ^ ^^^^ ^ _ ^2.23) 



32M 



AT 



In (|2.23|) , we used the Goldberger-Treiman relation q^^^^/Mm = Qa/ fn- Note that in (|2.19| ) 



the two structures proportional to 1 and 4 have opposite signs while they have the same sign 
in ( p:23D . 



One can easily check that the crossing term between (a) and (b) disappears. 
In the above (a) and (b), we have obtained the ir + N ^ N scattering contribution to 
the spectral function at q = in the form: 

p-+^-^(.;,T) = pJ^^''-''{q,T)+pl^'^-^{q,T) 
( A^C ■ 



32M, 



N , 



{5{uj - Mm) - 6{u + Mm))[{i + qI) ^ - (l - (?^) M^} . 

(2.24) 



We note that up to 0{T'^) the use of the nucleon mass and the pole residue of T = in 
( |2.24| ) is consistent with our treatment of the OPE side because of the presence of the factor 
C in (|2.24|) . (From a physical ground, one might wish to replace them by those at T 7^ 0. 



But these two procedures cause only an O(T^) difference in the final result.) 

In (|2.24|) , the n + N ^ N scattering term eventually becomes an effective delta function 



at the pole position of T = 0. However, it contributes differently to the 1 and 4 structures 
of the correlators. Therefore its effect must be taken into account in the spectral function 
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when we make use of the Borel sum rule method. Otherwise, an erroneous mass shift of the 
nucleon would occur. 

(ii) vr+A^ — > A contribution: Using the lowest order piece of the chiral invariant 7c~N — A 
interaction lagrangian Ct^na ~ Q-nNA^t^l^-^^^'^ HDI; calculate the vertex: 



(7r(A;)|r^(0)|A(p)) 



N 



(2.25) 



where A^(p) is the Rarita-Shwinger spinor for A. The sum over spin for A gives the projec- 
tion operator ||20l 



spin 



1 



3Mi 



1 



3Ma 



{Pf^lu - Pul^) 



+ Ma).(2.26) 



In the soft pion limit, the 7r + — A scattering term appears at ~ Ma where the nucleon 
propagator in (|2.25|) becomes proportional to 1/{M\ — M^) and the vertex contribution is 
k^k^Vfj^u ~ k"^- Therefore the vr + — ^ A scattering contribution becomes O(T^) and we 
shall discard it. 

One can easily repeat the same steps as above (i) (a) (b) and (ii) for the other octet 
baryons. By the same reason as (ii), the effect of the transition to the decuplet baryons 
(tt + S ^ S*(1385), vr + S ^ H*(1530)) is O(T^). Even among octet baryons, (i)(b) type 
scattering between A and S is 0(T*^) because of their mass difference. 



2.3 Borel sum rule 

We are now ready to perform the Borel sum rule analysis for the octet baryons. We assumed 
that the finite-T medium is the dilute pion gas which has zero chemical potential, and thus 
the charge conjugation symmetry is preserved, i.e., there appears no splitting between the 
baryon and anti-baryon poles.0 Therefore the spectral function for the nucleon reads at 
g = 

p^(«,r) = Xl{T)5{u'-Ml{T)){4+M^{T))sign{u)+p-+''-''{u,T) 

+ (continuum by step function), (2.27) 

where p^+^^^ is defined in (p.24|) . Corresponding to ( p.2|) , we decompose the spectral 
function as p^{q) = Pi{q) +4P'ii.'i)- Then both pi and p2 satisfy the relation pj(— cu, q = Qi) = 
—pi{ijj,q = 0) (z = 1,2) as was obtained in the previous subsection. In the deep Euclidean 
region, the dispersion relation ( |1.6D can be written as 

nf (^2 = -Q\ q = 0,T) = ^,(Q^ T) = r dn^^^^TT^^- = 2) (2-28) 

^ This is in contrast to the system with a finite baryon number such as the nuclear matter. In order to 
organize sum rules for baryons in such a system, the dispersion relation needs some modification. See |pT| 
for the detail. 
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Applying the Borel transform 

= lim 



QVn = M2 : fixed ) 



{Q 



2\n 



n,{Q^,T) (2.29) 



to ( p.9|) and (|2.1CI| ) with the spectral function ( p.27|) , we get the following relation: 



c 



2a 1 - ^ = ~X%{T)Mr,{T)e-^W/^' - 



16 



e 



(0)/M2 



(2.30) 



(2.31) 



+ M% + -a^ = A^(T)e-^^-(^)/^' + ^^(0)(1 + ^/)C ^-a4(o)/m^^ ^2.32) 
3 16 



where 

a = —{27rY{uu), 

TT 

~XUT) = 2(2vr)%(T). 
From ( |2.31D and ( ^.32| ) one gets the expression for the nucleon mass: 

2a (l - I) + (A?,(0)(1 - g/)Mj,iO)C/l6) e-*^^/^' 



M6 + M26 + |a2 - (A^(0)(1 + g/X/w) e 



A/2(o)/M2 



(2.33) 
(2.34) 

(2.35) 
(2.36) 



In ( 2.36 ), we omitted the correction due to the continuum contribution for brevity. The 
formula which include the correction is given by the following replacement: 



M^/l- fl + A , 

\ \ 2MV 

M^{l-(l + |L)e-o/M^} 

M2 (l - e-^«/*^' 



(2.37) 



An important consequence from ( |2.36| ) is that Mn{T) is completely T-independent. In fact, 
if we replace A^(0) in the numerator of ( p.36|) by the one obtained from ( |2.31| ) by setting 
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T = 0, and replace A^(0) in the denominator of (|2.36|) by the one obtained from (|2.32|) , 
we can easily see that the T-dependence disappears from Miy{T). If we did not include the 
scattering term, the nucleon mass would behave as M^iT) = M^r (0)(1 ~ (/8) as can be 
seen from ( p.36|) . This means that the T-dependence of the condensates caused through the 
interaction with the thermal pions is completely compensated by the it — N scattering term 
and the pole position does not move at least to order O(T^). This is consistent with the 
statement of Leutwyler and Smilga based on the chiral lagrangian P]: If we calculate the self 
energy of the nucleon using the iiNN effective lagrangian ( p.21| ), we can easily check that the 
real part of the self energy is zero at the T = pole position. This is a direct consequence of 
the Adler's consistency condition required for the use of PCAC (in another word, ( |2.21| ) 
has a derivative coupling). 

The pole residue changes as X%(T) = A^(0)(1 — (1 + g'^)(/16) as is seen from ( p.31|) and 
( |2.32| ); the same result obtained by in the chiral lagrangian approach. 

From the above demonstration, it is clear that we have to take into account the new 
structure in the spectral function consistently with the T-dependence in the OPE side of the 
correlator. Otherwise the usual procedure in the sum rule leads to an artificial change of the 
resonance parameters. 



In [10|, the nucleon mass was calculated by a finite-T QCD sum rule method. The authors 
found a dropping nucleon mass even in the low temperature region as {uu)t decreased at 
finite-T. They did not take into account the tt + N ^ N scattering effect and assumed 
that the T-dependence of the four-quark condensate is the same as {{uu)tY. Although their 
calculation was not based on the pion gas approximation, the present consideration shows 
it is crucial to treat both the OPE side and the phenomenological side consistently. Correct 
treatment of the T-dependence of all the condensates is also required. 



3 Analysis of Correlators using PCAC 



In this section we shall examine the octet baryon correlators starting from the pion gas 
approximation without using OPE: 



n(g,T) ~ n(g,0) +i / d'^xe"'"' 



(27r)32A;0 



nsikyT){7T'^ik)\T{vix)m) W^ik)). (3.i; 



Applying the LSZ reduction formula for the pion, then using the PCAC relation 9'^yl^(x) 
f.,^m'^(f)°- (x) for the pion field (j)°'{x) and taking the soft pion limit, one arrives at 



n(g,T) ^ n(g,o) 



24 



+ {t{[t^{x')M4 [:F'M.m])) + (T([^5^(x°),r/(a;)] [^5^(0), r^(0)]))} 
+ ■-, (3.2) 

where denotes the terms associated with the axial charges carried by the octet baryons 
(terms with in the previous section), which becomes O(T^) or higher except at the pole 
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position of T = O.Q Utilizing the formula (|2.17|) together with 

3 



srab 



5 ' 



T" n 
5' '/ 



N 



we obtain the following expression for the nucleon correlator at T 7^ 0: 
n^(g,T) = 



(3.3) 



(3.4) 



The second term of the r.h.s. of ( |3.4| ) corresponds to the tt + N ^ N scattering term in 
the sum rule approach in the previous section and the first term of the r.h.s. of (|3.4|) is the 
modification of the residue of the nucleon current. The above analysis of the current simply 
tells us that the nucleon correlator at T 7^ can be written as a superposition of the same 
correlator at T = with T-dependent coefficients and there is no O(T^) shift of the pole 
position, which is consistent with the observation of [Q. (As was shown in |^ using ( |2.21| ), 
the contribution to the real part of the self-energy from + ■ ■■ in ( p.4|) becomes zero at the 
pole position but induces an O(r^) wave function renormalization.) 

: n^(g,0) + ---. The 
and ( glbl) . This is 

the reason we observed no mass shift in the Borel sum rule. We remind the readers once again 
that the factorization of the four-quark condensate in the medium level {{qTqY)rp 
is not justified. If we adopted this procedure in section 2, the T-dependence of the OPE side 
of the correlator would be different from (|3.4]) . 

It is easy to extend the above analysis to other octet baryon correlators. For this purpose 
we need the following commutators: 



Equation (U) reads Hf (g,T) = Hf (g, 0)(l-C/8) + --- and n^(g,r) 
T-dependence of these two equations is the same as the OPE given in 



■a 



5 ' '/ 




21 



^ r + IT r + 



12 



V2\ 



r" 

2/11 



,(3.5) 
(3.6) 
(3.7) 



5 ' 



^ 5 5 



5 ' 



n 



where we introduced a new current 77^' defined by 



3 = 



7] 



(3.8) 
(3.9) 

(3.10) 
(3.11) 



■^In the sum rule, we needed to integrate over the spectral functions, which is why the 0{T^) contribution 
appeared. 
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in ( p.6|) and ( p.8|) . By using these relations in (|3.2| ), we obtain the following relations among 
the correlators: 

U\q,T) = n^(g,0) + ^n^^'(g,0)-^75n^(9,0)75 + ---, (3.12) 



n^(g,T) = (i-^)n^(g,o)-^75n^'(g,o)75 + ---, (s.is) 



n^(g,T) = |'l-l'|n=(g,0)-^75n-(g,0)75 + ---, (3.14) 



n^^'(g,0) = z / d'xe^'^^{T{r^ti^)r^^'iO) + ^ai^HiO)))- (3-15) 



with 



The current rj^' is anti-symmetric under the exchange between u and d quarks and thus rj^' 
(or "fst]^ ) should have some overlap with isoscalar strangeness=-l baryons such as A(1115), 
A*(1405) {J^ = l/2~) as well as the vr — S continuum contribution. As is shown in |16|, there 
are 5 independent interpolating fields without derivatives for A. One can see by the Fierz 
rearrangement between s and u (or d) that r]^ consists of r]^ and those others. We tried to 
identify its structure by the vacuum QSR including the operators up to dimension-6, but it 
does not seem to have a dominant pole contribution. In any case, ( 3.13|) tells us that the 



finite-T S-correlator can be written as the modification of the residue and the mixing with 
A' correlator. In principle, the second term of ( p.l2|) also describes the modification of the 
residue. As for S, the situation is completely parallel with the nucleon. The T-dependence 
of the OPE expressions ( |A.15| ) and ( |A.16| ) is consistent with ( p.l4| ) as they should. 



4 Summary and Outlook 

In this paper we have presented an explicit demonstration that the O(T^) dependence of 
the condensates which appear in the OPE of the octet baryon correlators is totally absorbed 
by the scattering terms n + B ^ B' and the modifications of the pole residues, in the 
framework of the QCD sum rules. This result is consistent with the statement by Leutwyler 
and Smilga P| . The result stems from the fact that the baryon correlators in the thermal pion 
gas can be written as a superposition of the correlators of T = with T-dependent coefficients 
up to O(T^). The procedure for achieving consistency with this relation is somewhat intricate 
in the QSR. Therefore one has to pay particular attention to the consistency between the 
assumption made to estimate the T-dependence of the condensates and the new structure 
appearing in the phenomenological spectral function. 

In the chiral lagrangian approach, the O(T^) mass shift of baryons is caused by tadpole 
interactions such as rrisKKBB / {K is the kaon field) In the m„ = = 0, ^ 
limit studied in this work, the kaon or r] field always accompanies the tadpole contribution 
as in the case of the above interaction. Thus without including kaons or rj in the heat bath, 
baryons do not receive a 0{T'^) mass shift. However, the presence of those massive excitations 
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is suppressed as ~ e For example, the above term in the effective lagrangian causes 

the mass shift of the order of m,TV(24/2)5i(my^/T) ~ 1 MeV at T = 150 MeV. If we 
included the kaons and ?7's in the heat bath together with the nonzero strange-quark mass, 
the T-dependence of the OPE for the hyperons shown in ([A.11|) - ([A.16|) would be different 
from ( p.l2| )-( p.l4| ). This would lead to the O(T^) mass shift of the octet baryons, although 
it should be tiny because of Bi{mK/T) (~ 0.06 at T = 150 MeV). 

To go beyond 0{T'^), one needs more information on the pion structure functions (twist- 
4 effects), O(T^) or higher T-dependence of the scalar condensates and the analysis of the 
structure in (|2.2|) as well as more involved treatment for the phenomenological side (such 
as octet — >• decuplet transitions). These issues are beyond the scope of the present study. I 
hope the lesson we learned through the demonstration in this work will be useful for more 
advanced studies on these higher order effects. 
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Appendix 

Here we will present the nonfactorized form of the four-quark operators which appear 
in the OPE for the octet baryon currents. Nonscalar four-quark operators are associated 
with the twist-4 contribution in the pion structure function. We ignore them since they 
become O(T^) effect in the QCD sum rules at T 7^ 0. We decompose the contribution of 
the four-quark condensate into the two pieces corresponding to the two structures of the 
correlators: 

n4":ta;;(^) = nr'"- + n^'''"-^- (a.i) 

(i)iV: 



n^(?) = -if 



0, 






-1 










){ 






'5 






,2 



(A.2) 



3 
2* 



(A.3) 



where A" is the SU(3) color matrix and we explicitly kept the Nc{= 3) dependence. As is 
seen from ( |A.3| ), the four-quark operators always appear in the form of 

1 - ^) {qTqq'Tq' + •■} - ^{gFA^TA'^g' + -}■ 



We will henceforth use the abbreviation — |{with A"} to denote the second contribution, 
fii) A: 



3g2 



1 — —— \ ^4:{uudd) — {ua ^jjuda fj_yd) — A{uddu) + {ua ^^dda ^yu)^ 



- {with A'^} 



+ 



3g2 



1 - l^{uuss) + (ddss) - ^{wy^usj^s) - ^{d'y^dsj^s) 



-^(«7m75MS7/,75s) - ^(rf7/,75rfs7/,75s) 
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1 



{with A"^} 



(A.4) 



-2 
3^ 



3 5 - 3 - 



--{with A"} 

+ ( 1 - i{uuss) - {UJ5USJ5S) + ^(m7m«s7m^) - J('"7m75'"S7^75s) 



+{ddss) - {dj^dsj^s) + ]^{d'y^ds^^s) - ]^{d'y^-f^ds-f^-f^s) 



{with A"} 

1 \ r5 



2\}' W] (2^"^"^^^''^^ ^ ^ W/.75(^(^7;.75ii)| + ^ {with A"} 



(A.5) 



--{with A"} 

+ ( 1 - ^) {(m7mM57^5) - (m7m75M57^75s)} - ^ {with A"} 



(A.6) 



-1 



1 - 



1 \ r5 



+ (1- 



NJ 12 

1 \ ri 



NJ 12 



3 11 
{u-f^us-f^s) + -('ii7M75^iS7M75s)| - 2 {^^^ ''^"} 

(M') - ^(fe^)') - ^((^7^75^^)') + \{{ui^,uf)^ 



--{with A"} 



(A.7) 



1 - |-3(mmss) + {ua^^yusa^^s) + ^ (Ma/,,,75MScr/,,,75s) 



--{with A"} 



(A.8) 
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-1 

+ (1 



{with A"} 



~ {with A"} 



(A.9) 



Factorizing these four-quark operators into the square of the chiral order parameter, one can 
easily get the form of (|]|), (|Il0| ) and the following (^ll|)-(|Al|) at T = 0. To get the T- 
dependence of the four-quark condensates, we need to calculate the double commutators of 
(|A.2|) - (|A.9|) with the isovector axial charge. The calculation is tedious but straightforward. 
The following formulas are useful in carrying out the calculation (g = (u, d)): 



[•^5,^75?] 



:FtAl^r'q 

^l, qcT/iuT^q 

^5 ' Q(^i^iyl5r''q 
H.Ql.r'q 
^5,1Wr'q 



-975^ g, 

-qaf^^-f^T^q, 
-qc^t^uT'^q, 

[^5, Wsg] =0, 

-S^'qi^q, 
-6^'qq, 

-S^'^qaf.^q, 



(A.IO) 

After calculating the double commutators, we end up with other four-quark operators. Ap- 
plying the factorization to these four-quark operators, we eventually got (|2.9| ), ( p.lO|) and 
the following (^AlD-(^l|): 

^ U{uu) (l-i)- (ss)] qMQ') + ^.qMQ') 



nf(g,T) 



127r2 



967r^ 



n^(g,r) 



Am 



'uu) 



3g2 
-1 



1 H -{uu){ss) 1 



647r^ 



127r2 



A{uu) (1 - M -3(ss) ) ln(g2 



(A.ll) 



16 



nf(g,r) 



47]-2 



(ss)g2ln(g2 



327r 



3g2 



647r4 
1 ,a 



47r2 



327r2 ' TT 



_ r 

3g2 I 8, 



(A.13) 



(A.14) 



nf(g,r) 



47r2 



-1 



647r4 



1 



327r2 



TC 



2(£s)2 
3g2 



(A.15) 



(A.16) 



The T-dependence of (|A.15| ) and ( |A.16| ) is the same as ( p.l4|) . We also note the T-dependence 
of the four-quark operators is different in the different channels. 
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Figure Captions 

1 n + N ^ N scattering term in which rj^ couples to tt directly. 

2 TT + ^ scattering term in which rj^ couples to the nucleon that interacts with 

TT. 
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